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1 Abstract 

O 

It is shown that the Brillouin zone integral for the interstitial KKR- 
Green function can be evaluated accurately by taking proper care of 

H : 

the free-electron singularities in the integrand. The proposed method 
combines two recently developed methods, a supermatrix method and 
a subtraction method. This combination appears to provide a major 
improvement compared with an earlier proposal based on the subtraction 
method only. By this the barrier preventing the study of important 
interstitial-like defects, such as an electromigrating atom halfway along 
its jump path, can be considered as being razed. 
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I. INTRODUCTION 

The Korringa-Kohn-Rostoker (KKR) Green function methodEHi has proven to be a pow- 
erful and elegant approach to calculate the electronic structure of defects in metals@. The 
expressions to be evaluated are exact results of multiple scattering theory and the method 
has been applied successfully in calculating effects of charge transfer and lattice distortion, 
both on the electronic structure^ and on physical quantities like the Dingle temperature^ 
and the effective valency of migrating atomsi. Until recently its elegance was believedi'i 
to arise only after applying the muffin-tin approximation to the atomic potentials. This 
amounts to an exact description of the electronic structure within non-overlapping spheres 
only, being embedded in an average constant potential, called the muffin-tin zero. Fortu- 
nately it could be proven0~0 that the muffin-tin approximation is not necessary, opening 
the possibility to do full-potential calculations in the framework of multiple scattering theory 
as well. 

There remains one drawback of multiple scattering theory, and this will be the subject 
of the present paper. The muffin-tin zero or free-space reference system still appears in 
the expressions to be evaluated, since free-electron poles are present in the integrand of the 
KKR-Green functions. These plaguing singularities have to be handled with care. As far as 
substitutional defects are concerned this problem was solved recentlyliiJii by implementing 
a supermatrix method. However, for interstitial defects, such as hydrogen in metals and an 
electromigrating atom halfway along its jump path, the problem has not yet been solved. In 
this paper we want to present a solution. It appears that the supermatrix method formulated 
for substitutional defects can be extended to the interstitial problem. Supplemented with a 
subtraction method the expressions become manageable and evaluable to a high degree of 
accuracy. 

The paper is organized as follows. In section || the KKR-Green function matrices of 
interest are defined and the different existing computational approaches are reviewed briefly. 
In section [III] the supermatrix method will be presented. In section [TV] the subtraction 
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procedure will be described. In section |V] the subtraction method will be tested. The paper 
ends with some conclusions and prospects. 



II. RELEVANT MATRICES AND DIFFERENT APPROACHES 

In a calculation of the electronic structure of dilute alloys by use of the KKR-Green 
function technique two matrices show upEl, one for defects at substitutional sites 

T"' = -^—[ ^ke^n'M-^k) (1) 

ilBZ JBZ 

and one for an interstitial defect 



7 V ' = [ rf 3 fce ik - R ^'6 p (k)M- 1 (k)6 p ' T (-k). (2) 



The integrals run over the Brillouin zone (BZ) with volume Qbz- A lattice vector Rj is 
denoted by a label j and Rj-y stands for the difference vector between the sites j and j'. 
Arbitrary sites, including non-lattice sites, are indicated by a label p. All matrices are a 
function of the energy E and carry (suppressed) angular momentum labels L, which stands 
for (Zm). Both integrands contain the inverse of the KKR matrix M(k) given by 

M(k) = r 1 - 6(k) (3) 

in which the matrix t expresses the scattering properties of a host atom, and for spherical 
scatterers is equal to — sin S e zS , 5 being scattering phase shifts to be labeled by the angular 
momentum label /. The matrix 6(k) follows from the matrix ^(k), defined by 

feP(k) = ^B ra e-*- R ™, (4) 

j 

after taking for p a lattice site label. The matrix ^(k) is the Fourier transform of the free 
space propagation matrix element B v \ given by 

B% = Aixi^'-^^'CLv^ht,,^), (5) 
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in which Clul" are Gaunt coefficients and h~£(r) = '(/cr)li(r) . Real spherical harmonics 
Yi{f) are used, hf are spherical Hankel functions and k = yE. The electronic structure of 
the metallic host follows from the condition 

det M(k) = 0, (6) 

which is a basic result of multiple scattering theory. 

The free electron singularities are hidden in the matrix 6 p (k). They are readily made 
explicit by writing down its reciprocal space representation 

^' (k) = ¥ e (k + K.)»-^ (7) 

in which 

F l( k) = S-y i( t)?M, (8) 

V 7TK M K3; ) 

In this expression the free electron poles at the energies E = (k + K n ) 2 are clearly present, 
K n denoting a reciprocal lattice vector. The latter equality in fact defines the so-called 
free electron sphere. Note, that for p being a lattice site the exponential factor reduces to 
unity and that only for that case the second term contributes. The ji are spherical Bessel 
functions. As usual in Green function treatments, the energy carries an infinitesimally 
positive imaginary part, which is indicated by E + . At this point we want to remark, that 
the formalism discussed in this paper is currently^ applied at complex energies as well. For 
that slight changes in the notation are required. However, in all calculations one has to 
approach the real-energy axis somewhere, so that the pole problem shows up anyhow. It 
is clear that the matrix M(k) also contains the free-electron singularities. At the k points 
defined by Eq. fl6|), which pertain to the electronic structure, its determinant value is zero, 
while at the free electron sphere it is singular. 

Looking at the integrand of the matrix Q pp ' it is seen, that it is full of singularities. Free 
electron singularities are present in the matrices ^(k) and 6 P ' T (— k). These are partially 
cancelled by those in the matrix M(k), but this matrix itself introduces poles corresponding 
to the electronic structure of the metal regarding the condition Eq. @. 
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In order to state the problem in actual calculations clearly we expand a little upon it. A 
popular way to evaluate the matrices and Q pp ' is using a subdivision of the Brillouin zone 
into microvolumes, cubes or tetrahedrons, going back to Gilat and Raubenheimer0, Jepsen 
and Andersen0 and Lehmann and TautH The required matrix inversion is performed by 
using its eigenvalues A and eigenvectors V, 



Mil = (VA-Vt) = £ V Lq \-Vt, q (9) 

q A q 

iggested by Lasseter and SovenEl and elaborated by Coleridge et al.0. The matrix 



as su£ 

T-"' is the simpler one and therefore has been calculated most intensively. This explains 
the availability of quite exhaustive studies of substitutional alloys and the relative lack of 
results for interstitial alloys§. 

As an introduction we concentrate on the different evaluation methods for T J ' J ' . These 
methods can be distinguished by tracing the influence of the infinitesimally positive imagi- 
nary part added to the energy E. Working out this influence explicitly a real and imaginary 
part of the integral in comes out according to the well-known equality 

_ = V- - niS(x) (10) 

X X 

The imaginary delta function part leads to a reduction of the Brillouin zone integral to 
an integral over a constant energy surface. If one is interested in properties at the Fermi 
energy Ep that integral runs over the Fermi surface. The tetrahedron grid reduces to a 
grid of triangles over the constant energy surface. The real principal value part remains. In 
practice it is evaluated along two different lines. The most straightforward line is to evaluate 
the corresponding Brillouin zone integral explicitly. This is achievable if one needs T JJ at 
the Fermi energy only, because quite a fine grid of k points is required. Coleridge et al.0 
proposed to use a weighted grid, having a finer subdivision of tetrahedra at the singular 
surfaces, the Fermi surface and the free electron sphere, and that is the way it is applied^. 
For electronic structure calculations one needs the matrix T- 73 ' at all energies, starting at 
the bottom of the band and going upwards. To that end another line of evaluation^ of the 



real part has been developed, taking advantage of the necessity to account for a considerable 
energy interval by employing a Kramers-Kronig relation. This relation expresses the real 
part of T JJ in terms of its imaginary part in the form of an integral over the energy. In 
this way the Brillouin zone integration is avoided at the expense of the necessity to evaluate 
the imaginary part, the constant-energy surface integral, up to relatively high energies. At 
the latter point some approximation has to be made, in choosing an upper-bound cut-off 
energy. It is worth while to note, that in practice in this method the free-electron problem 
does not enter in the integrals, neither in the evaluation of the real part, nor in integrating 
over the constant-energy surface. The integrand of the imaginary part simply is a product of 
unperturbed metallic-host wave function coefhcients§!. The free electron singularities only 
complicate the finding of the metallic bands, using Eq. (^|), and the corresponding wave 
function coefficients. 

A relatively recent developmentEl'0 is to evaluate the full Brillouin zone integral, sub- 
stituting E + = E + i8 and using a small value for 5. The grid of k points to be used is not 
as fine as required for the principal value part, because E + now is a complex number, while 
a value of 5 = 0.01E is already small enough to meet both the 5 — > limit and accuracy 
requirements. 

All this applies as far as the pole structure of M _1 (k) is concerned, corresponding to 
the metallic electronic structure. Now we turn to the singularities in M(k) due to the 
free-electron poles. These singularities do not look too serious, because M _1 (k) approaches 
zero at these points. This is so indeed, if one follows the so-called double linear method of 
Coleridge et al.0. In that method the full matrix M _1 (k) is represented as a ratio of two 
functions n(k) and d(k), and the integral I(T) over a tetrahedron T is written as 

I(T) = jf d 3 k^ = £ Ki{d x , d 2 , d 3 , d 4 K (11) 

The third member follows analytically from the second member assuming a linear behaviour 
of the functions n(k) and d(k) inside the tetrahedron. The numbers rii and di are the 
values at the four vertices. The weights Ki are given by Oppeneer and LoddeiH. The 



function d(k) is supposed to be zero at the singularities due to the electronic structure. In 
practice one takes for d(k) the eigenvalue of M(k) which becomes zero, so the one that 
determines the electronic structure. Denoting this eigenvalue by Ao, the function n(k) in 
the numerator is equal to Ao J2 q VhqVuql -V This holds for the alkali and noble metals. If 
the Fermi surface is composed of more than one sheet more than one eigenvalue becomes 
zero, of course at different k points, and a product of the corresponding eigenvalues is used. 
However, in this approach the function n(k) is not as linear as one might wish, even if just 
one eigenvalue becomes zero. In addition to smoothly behaving eigenvalues it contains the 
eigenvalue representing the free-electron singularity. Since, due to the inverse, n(k) at such 
points becomes zero the method still works, provided a relatively dense (weighted) grid is 
used. 

In a later development^l'0Ei it was considered to apply Eq. (|TT|) not to the full matrix 
M _1 (k), but to each term separately in the representation (|9|) instead. The functions d q (k) 
in the denominator are the roots X q , and n q (k) = Vhq^t'q- ^ * s cl ear that this cannot work 
without modification, due to the terms corresponding to the highly nonlinear free-electron 
roots. The modification requires an innovative handling of the free-electron singularities. It 
appears that the concept of a supermatrix has to be introduced in the descriptionilJii. This 
supermatrix method, which allows for a considerable reduction of the number of k points to 
be use , will be discussed in section |TJ. 

III. SUPERMATRIX DECOMPOSITION OF THE INTEGRANDS 

For the sake of clarity first the formulation for the matrix T-"' will be summarizedlll'0, 
after which it is given for the matrix Q pp ' . In the supermatrix method the terms in the 
matrix 6(k), which make the matrix M(k) singular, are treated separately in a special way. 
Suppose that a total number of iV reciprocal lattice vectors contribute to the singularity 
in the energy range of interest. Then, glancing at Eqs. (^|) and (0), it is possible to write 
the matrix M(k) as a sum of a smooth part M°(k) and a part containing the N possible 
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singularities as follows 



N 



M LL ,(k) = M° LL ,(k)-Y: 



F L (k + K n )F?,(k + K n ) 



(12) 



A 



After defining a square diagonal matrix D with elements D n = (k + K n ) 2 — E + and a 
rectangular matrix F with elements Fi n = ^L(k + K„), this equation can be written in 
matrix form as 

M = M° - FD~ X F^. (13) 

This form suggests inversion using a supermatrix A defined by 

D Ft 
F M° 

according to the Sherman-Morrison- Woodbury formulai! 

D- 1 + D^F^M^FD- 1 -D^F^M- 1 

A' 1 = 

-M^FD- 1 M" 1 

Contrary to the original matrix M, the supermatrix A is regular everywhere in the Brillouin 
zone. The free electron poles in M appear in its supermatrix representation A as free electron 
zeros of the matrix D in the upper left corner of A. It is even so that det A = det D det M, 
which is clear from the following equality 



(14) 



(15) 



det A = det 



D Ft 




( 


1 






= det 








F M° 




{ 


FD- 1 1 





D Ft 

M° — FD _1 Ft 
det 1 det 1 det D det M, 



(16) 



by using Eq. ([13|). One simple consequence is, that the unpleasant original KKR condition 
Eq. (H) can be replaced by the regular one 



det A(k) = 0, 



(17) 



This difference between A and M is crucial regarding the accuracy of the integration. Ac- 
cording to Eq. ([15]) the wanted inverted matrix M~ l is given simply by the lower right 
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block of the supermatrix A' 1 . In the inversion of A, to be achieved in a way similar to 
Eq. (|9]), using its eigenvalues and eigenvectors, all eigenvalues behave smoothly. Applying 
the double linear method symbolized by Eq. ([□]) to each term separately in the sum over 
the inverse eigenvalues, a mesh of about 100 k points is sufficient in most self-consistent 
electronic-structure calculations. The eigenvector products n q (k) in the numerator even be- 
have so smoothly, that it appears to be sufficient to use the value of the functions n q (k) in 
the middle points of the tetrahedron only. By that the third member of Eq. ([TT]) becomes 
proportional to that value, while the sum over the weights Ki reduces to one simple weight 
expression. A convincing example of the power of the supermatrix method is given in Ref.0 
Now everything is ready to focus our attention to the much more singular integrand of 
the interstitial KKR-Green function matrix Q pp ' . Fortunately a similar decomposition of 
that integrand can be designed. In addition to the form ([13]) for the matrix M one needs 
the forms 



Fik) = b m + F p D- l F\ 



(18) 



and 

b pT (-k) = b r0 + FD~ l F p] . (19) 
for the two b matrices in the integrand in Eq. (H). The matrix F p is defined by 

Fin = e iK "- R ^Ln (20) 



Using the form (pl|) for the supermatrix A 1 it is readily seen that the supermatrix P, 
defined by the following product of three supermatrices 



(21) 



after multiplication gets a lower right block, which is precisely the matrix product in the 
integrand of Eq. @). The left and right supermatrices in Eq. fl2~T|) do not behave smoothly 





D Ft 




D Ft 


-l 


D 


p = 


b l0 + F p D~ l F^ 




F M° 




F b r0 + FD- l F p '^ 
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yet, but after some rewritings, to be given in appendix the supermatrix P obtains the 
form 



P 



D 

~F p D- l F p '^ 










D Ft 


-l 





+ 






F M° 




b r0 



(22) 



The second term contains the inverse of the supermatrix A, now multiplied from the left 
and right by a smooth matrix. Since A" 1 can be obtained in a smooth way, the second term 
is easily evaluable. Only the free electron poles in the lower right block of the first term 
still require special treatment. It can be shown, that a subtraction procedure presented 
recentlyB allows for a quick and accurate evaluation of that term as well. This will be the 
subject of the next section. 



IV. THE SUBTRACTION IDEA 

Originally^ the subtraction method was designed and applied in handling the free elec- 
tron singularities in the integrand of the matrix Q pp ' according to Coleridge^, using a 
weighted distribution of tetrahedra in the Brillouin zone. For the sake of clarity that ap- 
proach will be summarized first. 

The idea of the subtraction method is to subtract a function /(k) from the integrand, 
which is chosen such that the integrand gets free of the poles, while the integral over the 
function /(k) can be evaluated analytically and is added later on. Although it can be seen 
by inspecting the behaviour of the integrand of Eq. (0) at the free electron poles that the 
function 

V e i(k + K n ).R pp , F L (k + K w )F£,(k + K w ) 

\ (k + K n P-£+ • (23) 

cancels these poles in the integrand, an explicit proof will be given in appendix [FJ. This 
function cannot be integrated analytically. However, it is possible to manipulate the form 
such, that it retains its pole-cancelling property on the one hand and can be integrated 
on the other hand. First the functions -Fz,(k) given by Eq. (f|) can be used in a simplified 

10 



form by taking the limit k — > k, by which the Bessel function factor reduces to unity. This 
step may induce some oscillations around the free electron singularities, but at the singularity 
the limit holds exactly. Another step is the introduction of an Ewald-like convergence factor 
in order to improve the convergence of the summation over reciprocal lattice vectors. The 
final function /(k) obtains the form 

f LL ,(k) = 2QbZ V e'W'V Yl( ^ + K ") yz/ ( k + K ") c -(lk+K n |- K ) 2 /^ ^ 
iik n (k -f" K n J^ E^~ 

in which the Ewald parameter rj controls the convergence. Upon integration over the Bril- 
louin zone, the summation over reciprocal lattice vectors leads to an integral over all k space. 
The angular part of the resulting integral can be carried out 

[ dy^™'Y L {k)Y L ,{k) = ^Y. %V ' C LL>L»3v\kR vv ,)Y LII {R vvl ). (25) 

The integral over k still contains a free electron pole, but the principal value part can be 
evaluated using the equality 



k z — hi Jo \ k + k 2 I k — k 

-{k-nf/r, 



+-rjl(KRpp>)V / 2dk— . (26) 



The first integral on the right hand side is regular, while the second integral equals the 
readily available exponential integral Ei(E/rj), being defined by@ 

E 1 (x) = / —dt. (27) 

J X t 

In this way the evaluation of the matrix Q pp> connecting two non-lattice sites R p and 
Rp' has been made possible, although it has to be admitted that in the used weighted mesh 
sometimes over 4000 k points are required. 

Now we return to the supermatrix P, Eq. (|22|). From Eq. (pD| ) one readily sees that 
the matrix F P D~ 1 F P '^ in the lower right block of the supermatrix in the right hand side of 
Eq. ( |2"2"D has precisely the form (|23|). Note that in Eq. ( P5| ) the exponential factor e lk ' R pp' 
in the integrand of Q pp ' , Eq. (§), has already been included, while the lower right block of 
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P represents the matrix product in Eq. (fj) only. So the only real difference between the 
matrix function —F p D~ l F p '^ and Eq. fl23|) pertains to the summation over reciprocal lattice 
vectors. In Eq. fl23| ) all of them are included, while in P only the iV poles-generating ones 
occur. From above it is clear that the function /(k) of Eq. ( p4|) has the same properties 
as expression Eq. Q2"5| ) as far as the free electron poles are concerned. It can be concluded 
that precisely that function can serve in evaluating the remaining problematic Brillouin 
zone integral. The function —F P D~ 1 F P '^ can be calculated straightforwardly, and after 
subtraction of /(k) given by Eq. (|2~4l) it becomes smooth. 

By this the barrier in evaluating the interstitial KKR-Green function Q pp ' can be con- 
sidered as being razed. In addition, the achievement of the supermatrix approach, in that 
a coarser mesh suffices compared with the original Coleridge approach, remains. The sub- 
traction method will be tested below. 



V. TEST CALCULATIONS. 

In this section the accuracy of two integrals will be tested, both of which suffer from the 
presence of free electron singularities in the integrand. The integrations will be carried out 
using the Coleridge approach^ symbolized by Eq. (|iTl). The first integral is given by the 
left-hand side of the following algebraic equality 



1 



b LL ,(k)^ R v>d 3 k = BH,. (28) 



Q B z Jbz 

The B matrix in the right-hand side, for different site labels defined by Eq. (|5|), is evaluated 
routinely up to any desired accuracy. If j' = j this matrix is defined to be zero, in accordance 
with the exact result for the integral on the left-hand side. The b matrix in the integrand, 
being given by the matrix b p of Eq. (|J) for p = j, is seen to be singular on the free electron 
sphere by inspecting the alternative representation, Eq. (|7|). Comparing with the matrix 
function F P D~~ 1 F P '^ in Eq. fl2"%D, clearly the singular part of the matrix b coincides with it 
if p and p' refer to lattice sites labeled by j and j'. So the integral of the b matrix can be 
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considered as a close test of the subtraction procedure proposed. Furthermore, the choice is 
quite natural, because this matrix is readily available in our computer codes. 

The actual tests have been done for copper at the Fermi energy. Using a muffin tin radius 
of 0.65a/2 and a lattice constant a of 6.831 Bohr, the phase shifts So to £3 are -0.1506388, 
0.0563578, -0.1491734, and 0.0010149 respectively and E F = 0.634 Ry. In the first test only 
the lattice constant and the Fermi energy enter. 

In Table H results obtained without and with subtraction are compared with the exact 
results according to the right-hand side of Eq. (|2q) . Only some representative matrix 
elements are shown, indicated by (j'/m)-labels in the first six columns. The seventh to tenth 
columns are obtained by numerical integration, the last column gives the exact values. The 
first row at the top indicates the way in which the integration is performed, and the second 
and third row specify the mesh and the denominator function d(h) in Eq. ([TT|) respectively. 
So, the seventh column follows from straight integration. Application of subtraction of the 
function /(k), defined in Eq. (0), is indicated explicitly. As for the second row at the 
top, straightforward integration requires a weighted mesh, which, in the present example, 
is denser near the free electron sphere, and therefore is indicated by FES. After subtraction 
the smooth integrand allows for a homogeneous grid, indicated by 0. Nevertheless the 
smoothness is tested by doing the same calculation with the weighted grid. The number 
of k points is given in the last row of the table. The singular surface is the free electron 
sphere and the denominator function d(h) required for the straightforward integration is a 
product of factors (k + K n ) 2 — E, as many of them that vanish somewhere in the Brillouin 
zone. Interestingly, this latter product is precisely equal to detD, the determinant value of 
the matrix D in the upper left block of the supermatrix A, introduced in Eq. (|i~4"|). This 
is indicated in the third row at the top. After subtraction d(k) can be chosen freely, and 
a constant, indicated by 1, is used. As an implicit test of the linearity of detD, also after 
subtraction a calculation is done using that function as a denominator, given in column 
10. Of course, in that case the weighted mesh FES is required again. Comparing the 7th 
and 8th columns with the last column it is seen that application of subtraction leads to a 
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major improvement. Furthermore, the better result is obtained with a much coarser mesh. 
Applying the finer mesh after subtraction, as shown in the 9th column, does not lead to 
significant changes. In addition it can be concluded from column 10 that the linearity of 
det D is sufficient. 

Another test can be derived from the matrix Q pp ' , see Eq. (Q) , if the arbitrary site labels 
are replaced by lattice site labels. Then, because V(k) = V T {— k) = b(k), the following 
algebraic identity 

Qii' = -t^Sjf - B"' + rY^'r 1 (29) 

can easily be derived, using the relation (0) between M(k) and b(k). Both matrices Q and 
T can be evaluated numerically only. However, while the integrand of Q^' contains free- 
electron singularities, in the integrand of T-"' they merely appear as a nonlinearity in the 
eigenvalues of the matrix M(k), as described in section |I[ In the tests is evaluated only 
according to the original Coleridge approach^, including a subtraction procedure as well, in 
the way the present authors have described recently^. The matrix T-"' will be calculated 
by the supermatrix method also, for which less k points are expected to be required. So 
features of both the subtraction and supermatrix methods will be illustrated by Eq. (|2"9"|). 

In Table |1| results for Q^' according to the left-hand side of Eq. (|29|) , before and after 
applying the subtraction procedure, are compared with results according to the right-hand 
side. This table is built up the same way as Table |. First we elucidate the denominator 
functions d{k) of Eq. flTT|), as specified in the third row at the top. As far as the poles 
corresponding to the electronic structure are concerned, the product of the vanishing eigen- 
values of M(k) can be taken, as described in section |I[ Det M(k) might work as well, being 
equal to the product of all eigenvalues. But the latter choice introduces a problem because 
of its singular (and therefore nonlinear) behaviour at the free electron surface. Therefore, 
det M(k) has to be multiplied by det D. By that additional multiplication, regarding Eq. 
(PS|), the denominator function d(k) becomes equal to det A, which appears in the modified 
KKR condition Eq. (|17])- The free-electron singularities in the integrand of require mul- 
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tiplication by an additional factor of detD. So for straightforward integration, see column 
7, det A det D has to be used. In that case the (weighted) mesh must be taken to be dense 
near both the Fermi surface and the free electron sphere. This mesh is indicated by FS/FES. 
After subtraction the free-electron singular surface is absent. This allows for the less dense 
mesh FS and for the simpler denominator function det A, column 8. To see the influence of 
the density of the mesh the same calculation is done with the denser FS/FES mesh, column 
9. In addition, also in this complicated case, the effect of the (extra) nonlinearity of the de- 
nominator is investigated, by doing the same calculation with the denominator det A detD, 
as shown in column 10. Columns 11 to 13 give results according to the right-hand side of Eq. 
(p9p, the last two columns showing supermatrix results for two different meshes. Contrary to 
the weighted meshes of the Coleridge method, the meshes used in the supermatrix method 
are always uniform. 

Again the results improve largely by applying subtraction. In fact, the straightforwardly 
evaluated integrals can be considered as really bad. An obvious source of inaccuracy is the 
behaviour of the denominator function det A det D in regions where a tetrahedron is cut by 
the Fermi surface as well as the free electron sphere. Then the product of det A and det D 
does not change sign and the denominator does not become aware of crossing a singular 
surface at all. This even may happen for a noble metal. The single sheet Fermi surface 
intersects the free electron sphere in the neck region. Comparing columns 8 and 9 one sees 
that the denser mesh near the free electron surface yields not too large but yet non-negligible 
modifications. The addition of the extra factor det D is not as harmless as in the integration 
of b(k), as can be seen from column 10. This is just another illustration of the source of 
inaccuracy indicated above. So it can be stated that it is important to use a denser mesh 
at the free electron sphere, even if the integrand is not singular there. This certainly is a 
demerit of the Coleridge method. Nevertheless, it is satisfactory to see (columns 11 to 13) 
that, in evaluating the T matrix, the Coleridge method leads to the same results as the 
supermatrix method. In addition the table confirms explicitly that the subtraction method, 
proposed recently^ for determining the full Q matrix, is reliable. 
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Finally, the table shows implicitly, that the supermatrix method for determining the 
Q matrix is more efficient than the subtraction methodlll. Less than 1000 k points are 
sufficient. The number of 891 in the last column is a good indication for integrating the 
second term in Eq. which contains the same supermatrix A^ 1 as the supermatrix 

integrand of the matrix T. The subtraction procedure proposed for the first term, having 
free electron poles only, is tested by Table |. The number of 640 k points in the 8th column 
suffices for that term. 



VI. CONCLUSIONS AND PROSPECTS 

The supermatrix method, initially proposedQ with the aim of a fast and accurate evalu- 
ation of the KKR-Green function appearing in calculations of the electronic structure 
of substitutional alloys, has been extended to the interstitial KKR-Green function Q vp ' . A 
subtraction procedure is shown to resolve the remaining problem of integrating a function 
with free electron singularities only. In evaluating the corresponding Brillouin zone integrals 
rather coarse grids of less than 1000 k points can be used. Applications to calculations of 
transport properties such as the effective valence of electromigrating atoms are in progress^. 
Study of the electronic structure of the largely unexplored interstitial defects will be the sub- 
ject of future investigations. 
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APPENDIX A: 



We want to give a derivation of the final form Eq. ( p^) of the supermatrix P from its 
definition Eq. (|21|). Using the matrix A defined by Eq. ( |14D we first write the matrix P a 
little more compact 



D Ft 

b 10 + F p D~ l F^ 



A' 1 



D 

F b r0 + FD^FP'i 



(Al) 



Subsequently the supermatrix left of the supermatrix A" 1 is written in the form X l A + Y l 
and the supermatrix to the right is written similarly as AX r + Y r . In principal, the choice 
of the four matrices X 1 , Y l , X r and Y r is arbitrary, but a suitable one is 



X 1 



1 













FPD- 1 







(A2) 





1 







x r = 




Y r - 
1 1 











6 r0 



(A3) 



The elaboration of Eq. (|A1|) is now straightforward and ends up with the expression 





D 







A' 1 


-FP't 


p = 




+ 








-F p D- l F p '^ 








b r0 



(A4) 



Regarding the definition of the supermatrix A, Eq. (|14D, this is precisely the required 
form Eq. ( p2|) given in the main text. 



APPENDIX B: 

We want to show that the function (^) indeed cancels the free electron poles in the 
integrand of the interstitial KKR-Green function (Q). Apart from the trivial exponential 
factor e ^k ' :R ^ , ^ , ', the integrand, being a matrix to be denoted as a, can be written in the 
following form 
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a = (b 10 + F P D' 1 F J( )(M - FD^F^iV + FD^F^), 
in which Eqs. (|l"3"D, ( |18D and (|19|) are used. First we define a matrix c 

c = (M° - F J D- 1 F t )- 1 (6 r0 + FD^F^), 
being equal to the product of two matrices in the triple product matrix a. So 

a = b lo c + FW^F^c, 
while Eq. ( |B2|) can be written in the form 

FD^F^c = M°c - b r0 - FD^F^. 



(Bl) 



(B2) 



(B3) 



(B4) 



Three observations can be made. First, from its definition it is clear that the matrix c 
is a regular one, because the poles in numerator and denominator cancel. Secondly, as a 
consequence, only the second term in Eq. ( |B3|) contains the free electron poles. Thirdly, the 
left hand side of Eq. (|B4|) already has the form of that second term. Regarding Eq. (|20|) 
it becomes equal to that second term by multiplying Eq. ( p4|) from the left with the factor 
e lKn ' Rp . By that the third term on the right hand side of Eq. (|B4|) obtains the required 
form —F P D~ 1 F P '^ . Substituting the rewritten form of Eq. ( |B4|) in Eq. ( |B3|) completes a 
first proof. The merit of the second proof to be given is, that in addition it exhibits a nice 
linkH between the subtraction and supermatrix methods. 

To that end an auxiliary rectangular matrix q is defined by the relation 



Ft c = -F p,t - Dq 

by which Eq. ([B4|) reduces to 

Fq = -M°c + b r0 . 
Now Eqs. (|B5|) and (|B6|) can be combined in a supermatrix form 



(B5) 



(B6) 



D Ft 




Q 






F M 




c 




b r0 



(B7) 
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It is seen that the supermatrix A, Eq. fll4|) , enters the formulation. 
The solution of Eq. ( |B7| ) 



Q 


= A' 1 








A' 1 




__pP't 










c 




b r0 




A- 1 

A 2\ 


A' 1 

^22 




b r0 



(B8) 



requires the subblocks A^ of the supermatrix A 1 , by definition given by 



D Ft 
F M 

We only need the two explicitly written equations 



A' 1 A' 1 

^ll -^12 



A~ l A- 1 

^21 ^22 



/ 
/ 



(B9) 



F)A[l + F*A£ = I and DA^ 1 + F ] 4 1 = 0. 



x 21 



[2 



l 22 



After substitution of the matrix c according to Eq. ( |B8| ) in Eq. 
interest a, 



(BIO) 

for the matrix of 
(Bll) 



a = (b m + F p D- l F ] )(-A 2 ^F p ' ] + A^b r0 ), 

and making use of the equalities (P10| ), one finds 

a = -b l0 A^F p '^ + b l0 A^b r0 + F P A^F P '^ - F p A^b r0 - F P D~ 1 F P '\ (B12) 

in which form once again the term containing the free electron poles is made explicit. 

Up to now the matrix a figures in the subtraction method only, standing for the singular 
matrix in the integrand of the interstitial KKR-Green function matrix Q pp ' . It is interesting 
to go one step further and to write Eq. ( |B12| ) in the matrix form 



a = [-F p b 10 ] A 



b r0 



F P D~ 1 F p,J< . 



(B13) 



In this form one recognizes the lower right block of the supermatrix P, Eq. (|A4|) , which 
figures in the supermatrix method. So the supermatrix and subtraction procedures appear 
to be intimately linked. 
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TABLES 

TABLE I. B^ L , calculated by integration of 6(k) in different ways. The last column contains 



its exact value. 
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straight 


J( k ) 




Jl k ) 


exact 














r Hjo 


n 
U 




TT'TT'Q 




3 


3 


t 


m 


V 


TTl' 


aet 


i 
1 


i 
1 


Ant n 

aet 




1 


-t 
1 


n 

u 








I) 


n noo 
0.023 


O.OUO 


n nnn 
U.UOU 


n nm 
U.0U1 


n 







1 





1 


U 
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3 


o 
Z 


U.zzo 


n ni ^7 
0.01 1 


U.UU6 
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0.040 


0.000 


0.U00 


-0.002 





2 


1 











n 


207 


1 99 


n 199 


905 


198 






1 





1 





0.164 


0.172 


U.172 


0.178 


0.172 






2 


1 


2 


l 


0.464 


0.470 


U.47U 


0.477 


0.469 






3 


2 


3 


2 


2.794 


2.649 


2.647 


2.709 
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1 





3 





0.457 


0.456 


U.455 


0.457 
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Number of k points 




2247 


640 


2247 


2247 
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TABLE II. The matrix Q^ L , calculated along different lines. 





G 


on' 

LL' 


labels 




straight 


/(k) 


/(k) 


/(k) 


r.h.s. 


r.h.s. 


r.h.s. 












FS/FES 


FS 


FS/FES 


FS/FES 


FS/FES 


super 


super 


3 


f I 


m 


/' 


m' 


det ,4 det £> 


det A 


det A 


det A det D 


det A 


matrix 


matrix 


1 


1 
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-0.149 


-0.155 
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1 
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0.008 
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0.160 
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2 


1 


2 


1 


0.351 
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2 


3 
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2 
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-0.151 
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2 


1 


2 
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3 


2 


3 
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-0.027 
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-0.079 
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Number of k points 


4123 
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4123 


4123 


4123 


146 


891 



23 



